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Abstract. Using the Paul Trap as a model, we point out that the same 
wave functions can be variously coherent or squeezed states, depending 
upon the system they are applied to. 

1 Introduction 

Elsewhere [[Tl-IQI? we have investigated time-dependent Schrodinger equations that are 
quadratic in position, q, and momentum, p = —idq-. 

{2idt - [i + g2{t)]p^ + go{t)\{qp + vq) 

+ g^{t)p-2h2{t)q^ -2h^{t)q-2ho{t)]^{q,t) = Q. (1) 



We use dimensionless variables {% = m = 1) and the (7's and /I's are time-dependent 
functions. This study yielded methods for obtaining the symmetries of these systems. It 
also provided explicit analytic relationships between various subclasses of these equations. 
From that we could obtain implicit and sometimes explicit analytic solutions for the 
number, coherent, and squeezed states, and also for the associated uncertainty relations. 

^Our discussion of the physical distinction between coherent and squeezed states is most appropriate 
as a vehicle for us to pay homage to the late Dan Walls. Since the time of their experimental discovery 
in the mid 1980s, his name has been closely associated with the physics of squeezed states. 



In this paper, we will apply the techniques so obtained to the Paul trap system. We 
will thereby be able to provide insight into what are the coherent and squeezed states of 
the Paul trap, and into what uncertainty relations are satisfied. 

In Section 2, we will set up the formalism for the particular subclass of Eq. (|l]) that 
contains the Paul trap as a special case. That is: all hi{t) = gi{t) = except /i2(^) = /(^)- 
We go on, in Section 3, to describe the coherent states, squeezed states, and uncertainty 
relations for this subclass. The physical Paul trap is described in Section 4. 

In Section 5, we discuss coherent states, squeezed states, and uncertainty relations 
for the Paul trap. This is the main thrust of this paper. Specifically, the reader will see 
how the same physical states can be viewed either as coherent states of the Paul trap or 
as squeezed states of the harmonic oscillator, and vice- versa. 

2 Time-Dependent Quadratic Systems 

We now concentrate on the aforementioned particular subclass of Eq. (|l|), the time- 
dependent, harmonic-oscillator, Schrodinger equation. We have thatQ 

{d,, + 2idt-2f{t)q^}^{q,t) = 0. (2) 

Here the function f{t) is a continuously differentiable and integrable function of t. From 
the general results for Eq. (|I|), we have that Eq. (0) admits a symmetry alge- 

bra that is a product of a Heisenberg-Weyl algebra, wl, with a sm(1, 1) algebra. The 
subalgebras wl and sm(1, 1) are the ones associated with coherent states and squeezing. 
The time- dependent lowering and raising (ladder) operators in are given by 



Ait) = I [m p - iit) q] , A\t) = -I [r (t) P - C{t) q] . (3) 
Here q and p = —idq are the ordinary position and momentum operators 



a + a"'' a — a"!" . 



It is straight forward to demonstrate that the time- dependent functions, ^ and ^* in 
Eq. are two linearly independent, complex solutions of the second-order differential 
equation! 

7 + 2/(t)7 = 0. (5) 

Therefore, the time-dependent functions ^ and ^* are determined by the particular Hamil- 
tonian. 



^The Hamiltonian is H = —\dqq + f{t)q^. 

^ The 'dot' over a function of t indicates ordinary differentiation with respect to t. 
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The Wronskian of these solutions is 

mc{t)-mcit) = -t. (6) 

On account of Wronskian (|^), the ladder operators A{t) and A^lt) satisfy the commuta- 
tion relation 

[A(t),At(t)] = /, (7) 
where J = 1 is the identity operator. 

3 Coherent/squeezed states and uncertainty relations 

3.1 Coherent states 

First, consider \a;t), the displacement-operator coherent states (DOCS). They are de- 
fined by 

\a;t) = DAia)\0;t) = exp \aA^{t) - a*A{t)] \0;t). (8) 

In the above, DA{<y) is the displacement operator and a is a complex constant. Similarly, 
the equivalent ladder-operator coherent states (LOGS) are defined as 

A{t)\a;t) = a\a;t). (9) 

In Eq. (^), the extremal state |0; t) is a member of a set of number states, {\n;t), n = 
0, 1, 2, ■ ■ ■}, which are eigenfunctions of the number operator A''{t)A{t) [Q, 

A\t)A{t)\n;t) =n\n;t). (10) 

An important conceptual point for these time- dependent systems is that the number 
states, \n;t), are in general NOT eigenstates of the Hamiltonian. That is, 

H\n;t) =idt\n;t) ^ (real constant) |n; t). (11) 

In particular, the extremal state |0;t) is technically NOT the ground state. 

3.2 Uncertainty relations and squeezed states 

The operators, A and A^ , are linear functionals of q and p. With them, generalized 
position and momentum operators can be defined: 

^ A + At A-A^ 

Q= ^ , P = 12 

V2 2^2 

Their associated commutation relation is 

[Q,P]=iO = i[A,A\ (13) 
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where O is an Hermitian operator.^ 

Now consider the uncertainty product 



UiQ,P) = {AQY{APY = [{Q^) - {Q)^][{P^) - (P)^] (14) 

= m - (Qmmip - (p)?)]- as) 

Applying the Schwartz inequahty to Eq. (pTSf ) we have 

(AQnAPy > l\^{0) + {{Q,P})\' (16) 

> \{0)' + l{{Q,P})' (17) 

> (18) 

{, } being the anticommutator, and 

Q^Q-{Q), P^p^(p), (19) 



Eq. ( p!?!) is the Schrddinger Uncertainty Relation g]. Equahty is satisfied by states 
\B,C) that are cohnear in Q and P: 

Q\B,C) = -i B P\B,C), (20) 

where 5 is a complex constant. Now going to wave-function notation, for what will be 
the minimum-uncertainty squeezed states (MUSS), Eq. (pO]) can be rewritten as 

[Q + tBP]^Ut) = C^ssit), (21) 
C = {Q) + tB{P). (22) 

In general, B and C are both complex. The solutions to Eq. (|21| ) are "squeezed states" 
for the system. 

B can be understood to be the complex squeeze factor by i) multiplying Eq. ( POj ) on 
the left first by Q and then by P and taking the expectation values, ii) then doing the 
same for the adjoint equation (multiplying on the right), and iii) finally using Eq. (|15D . 
This yields 



B - ^^-7^, (23) 
{QP) (AP)2' 

- (24) 

IPP - (25) 



^ Because of Eq. (|^), for the systems studied in this paper O ^ I. But for now we continue with the 
general O. 
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Thus, \B\ yields the relative uncertainties of Q and P. 

For the particular case B = 1, the MUSS are the minimum-uncertainty coherent 
states (MUCS). These satisfy equality of Eq. (|T^), the Heisenherg Uncertainty Relation. 
That these are coherent states is easily seen from Eq. (^T]) , which then reduces to \/2 
times the LOGS Eq. (1). 

To intuitively see this, consider Eq. (^) for the simple harmonic oscillator. The 
solutions are 



ilJss{,q) oc exp 



2 



B 



i{p) q 



(26) 



B is an extremely complicated functional 
su(l,l) squeeze operator:^ 



of the complex parameter A of the standard 



Sa{X{t)) = exp [|A(t) a^a^ - |A*(t) aa 



\{t) = r{t)e"^^'\ 



(27) 



The exact relationship depends on if one defines the displacement-operator squeezed 
states (DOSS) as 

\a,X-t) = DA{a)Sa{\)\Q) (28) 

or as Sa{X)DA{oi)\Q) . 

We will now describe and then apply this formalism to a well-known and important 
system, the Paul trap. The aim is to obtain a deeper insight into the "coherent" and 
"squeezed" states of this system. 



4 The Paul Trap 



The Paul trap is a dynamically stable environment for charged particles It has 

been of great use in areas from quantum optics to particle physics. 

It's main structure consists of two parts. The first is an annular ring-hyperboloid of 
revolution, whose symmetry is about the x — y plane at 2; = 0. The distance from the 
origin to the ring-focus of the hyperboloid is tq. The inner surface of this ring electrode 
is a time-dependent electrical equipotential surface. The second part of the structure 
consists of two end-caps. These are hyperboloids of revolution about the z axis. The 
distance from the origin to the two foci is usually do = tq/ The two end-cap surfaces 
are time-dependent equipotential surfaces with sign opposite to that of the ring. The 
electric field within this trap is a quadrupole field. With oscillatory potentials are applied, 
a charged particle can he dynamically stable. 



Paul gives a delightful mechanical analogy ]rT| . Think of a mechanical ball put at the 
center of a saddle surface. With no motion of the surface, it will fall off of the saddle. 
However, if the saddle surface is rotated with an appropriate frequency about the axis 



^ Here, we can allow A = A(i). But note that if, in contrast to the present case, the oprators defining 
the squeeze operator have time derivatives (e.g., the su(f ,1) operators M_ and M+ of Ref. |^), then 
this possibility leads to time-derivative complications and an inequivalent result. 
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normal to the surface at the inflection point, the particle will be stably confined. The 
particle is oscillatory about the origin in both the x and y directions. But it's oscillation 
in the z direction is restricted to be bounded from below by some zq > 0. 
The potential energy can be parametrized as 



where 



V{x, y, z, t) = \4(x, t) + Vy{y, t) + V,{z, t), (29) 



^ „2 _ In /-i^\^2 



V,ix,t) = + Vit) x' = -n.,it)x\ (30) 
2rg 2 



Vy{y,t) = +^^V{t)y'^^Qy{t)y', (31) 



V.iz,t) = -^Vit) z'^ln^it)z'. (32) 



r, 







In the above, 

V{t) = Vdc - Vac cosiu{t - to) (33) 

is the "dc" plus "ac time-dependent" electric potential that is applied between the ring 
and the end caps. These potentials can be used to solve the classical motion problem. 
The result is oscillatory Mathieu functions for the bound case The oscillatory motion 
goes both positive and negative in the x — y plane, but is constrained to be positive in 
the z direction. 

Exact solutions for the quantum case were first investigated in detail by Combes- 
cure 



121 . In general, work has concentrated on the z coordinate, but not entirely [| 



Elsewhere we will look at the symmetries, separations of variables, and the number 
and coherent state solutions of the three-dimensional Paul trap, in both Cartesian and 
cylindrical coordinates. 



5 Coherent physics of the Paul trap 

5.1 Coherent states and the classical motion 

With the background established in the previous two sections, we want to discuss the 
coherent /squeezed states of the Paul trap. We focus on the interesting z coordinate and, 
in particular, use as reference the lovely discussion in Schrade et al. (SMSG) [1^. Using 
Eq. (P^, the Hamiltonian is 

H = ^^t = -l^ + ln^{t)z^ (34) 
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The connection between the notation of Section |^ and that of SMSG [15] is 

q = z, 2f = Q,, ^{t) = j^e{t). (35) 

The ladder operators, A{t) and A^{t), are 

A{t) = -^[eit)p-eit)z] = l{a[eit)-^m]+a^[-eit)-^eit)]}, (36) 
A^t) = -^[e*{t)p-e*{t)z] = \{a^[e*{t) + ze*{t)]+a[-e*{t)+ze*{t)]}, (37) 



Recall that because of Eq. (|35D , e(t) is a complex solution to the differential equation 
(|^) and its complex conjugate, e*(t), is the other linearly independent solution. The 
Wronskian of e(t) and €*{t) is a constant: 

W{e, e*) = e{t)e*{t) - e{t)e*{t) = -2i. (38) 

Note from Eq. (|3^), that Newton's classical equation of motion is 

F = Ut) = -^^^^ = ^ V{t) z ^ -Q,{t)z,i{t), (39) 

the solutions being Mathieu functions. But Eq. ( ^5|) shows that Eq. ( ^9|) is exactly the 
form of Eq. (^ for e{t) and €*{t). This means that a combination of e(t) and e*(t), up 
to normalization, follows the classical-motion Mathieu function solutions for Zci{t) and 
Pci{t) = Zci{t). Therefore, we can conveniently and with foresight write these classical 
solutions in the forms most convenient for the quantum study: 



^ci{t) = '-{[e*{t)e,-e{t)e:]p,+ [e{t)i:-e*{t)io]zo}, (40) 
Pciit) = '-{[t{t)e,-e{t)e:]p, + [e{tX~e*{t)^x,} = Ut), (41) 



where Zo and Po are initial position and momentum and eo = ^{to)- 
Now using Eqs. (0), (j^), and (^, with Q ^ Z, 



m = l{W)-^*it)]p-m-i*it)]z} (42) 
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[a [{Im e{t)) + i{Im e{t))] - a ^[{Im e{t)) - i{Im e(t))]} , 



Pit) = ^{m+e*it)]p-m + e*{t)]z} 

= {a [{Re e{t)) - i{Re e{t))] - a t[(i?e e{t)) + i{Re 



(43) 
(44) 
(45) 



It follows that 



-Me* 



(t)e(t)-e(t)e*(t)] = z, 



(46) 



\m.m\- 2 

where we have used the Wronskian (p^). 

Using Eqs. (H), (0), and (pT]), one can calculate that the coherent-state wave func 
tions for the Paul Trap go as 



^cs{,z,t) ~ exp 



2 I <t) , 



{z - Zci{t)f + izpciit) 



(47) 



The Gaussian form of \l/ can be readily verified @, |^ by noting that i) 



and ii) is a positive real function of t. 
5.2 Uncertainty relations 



m 

e{t)e*{t) 



1 - M) 



(48) 
(49) 



However, with these wave functions, SMSG ||T5| found that the z — p Heisenberg Uncer- 
tainty Relation is not satisfied. Rather, the z — p Schrodinger Uncertainty Relation is 
satisfied |15|, |I6|: 



(50) 



where (j) is given in Eq. (^) .f| 

^ The function (Es) also arises in number-operator states, the relevant equation is A(t)\0;t) = 0, 
when A{t) of Eq. (|36|) is used. This is the 0-eigenvalue case of Eq. (||). 

^ The right hand side of Eq. ( |50| ) is equal to 1/4 for all t only if = 0; that is, for an harmonic 
oscillator or a driven oscillator llIBi . 
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We observe that not only is this correct but is to be expected. These wave func- 
tions were generated by the Z — P variables, and hence should satisfy the Heisenberg 
uncertainty relation for Z and P. Contrariwise, note that Eq. (|36| ) can be written as 



A{t) = ^ {a[e{t) - ie{t)] + a^[-e{t) - ie{t)]} = fia + ua) . 



(51) 



But with Eq. ( p8D one has that 

1 = |^|2 _ |;y|2. (52) 

That is, A and must be related to a and by a Holstein-Primakoff/Bogoliubov 
transformation [|l^ of the form 

S-^{X)aSa{X) = [cosh r]a + [e^^sinhr] at. 



S-\X)zSa{X) = ^{ 



cosh r + e * sinh r 



a + 



cosh r + e* sinh r 



(53) 
i^}. (54) 



However, there remains one further complication. The coefficient of a on the right 
of Eq. (|53D is real while that of Eq. (|5TD is complex. There is a phase offset. This 
is related to the fact that there are four parameters in Eq. ( pT]) and only two in Eq. 
(|53|). That there are only two parameters in Eq ( pSf ) follows from the fact that the entire 
squeezed-state transformation from {z,p) to {Z,P) also involves a displacement. (See 
Eq. (pSl).) The displacement supplies the other two parameters, as we now demonstrate. 



5.3 Squeezed states 

Combine, in the A representation, the DOCS definition of Eq. (^ with the LOCS 
definition of Eq. (|^): 

A{t)DAia)\0;t) = aDA{a)\0;t). (55) 

Next, insert / = S~^{X) Sa{X) in front of all the operators and multiply on the left by 
Sa{X) of Eq. (|27|) . Regrouping, we obtain the equation 

{SaAS-') (SaDA{a)S-') Sa{X) \0-t) = a (^„Da(«)^-') Sa{X) \0;t). (56) 



Given Eqs. ( ^if) and (|5^) and using both the BCH relation 



exp 



A- 



exp 



7+- 



exp 



73(a^a + i) 



exp 



7-- 



aa 



(57) 



7_ = — e ' tanhr, 7+ = e* tanhr, 73 = In cosh r. 



(58) 
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and also the theorem 
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it follows that 



exp[X] exp[y] exp[-X] = exp[e^Fe 



Sa{X)DA{a)S^\\) = exp \pa^ - /J^al = D,{(3) 



/J * * /O* * * 

p = av — a u, p = a V — au . 



In addition, we see that 



where the coefficients u and v are 

u = u cosh r — fie^^ sinh r, v = fi cosh r — i/e~*^ sinh r. 
Note that u and v are functions of t because fi and z/ are. Furthermore, we have 

I |2 I |2 I |2 I |2 1 

|f I — \u\ = \fl\ — |z/| = 1. 
Therefore, Eq. (^) becomes 



(59) 

(60) 
(61) 
(62) 

(63) 
(64) 



{va + ua^} Da{f3) Sa{X) |0; t) = {a(A, /5)} Da{f3) Sa{\) |0; t) = {a(A, /5)} |a(A,/5),t). 

(65) 

This satisfies both the displacement-operator definition of squeezed states (to the right 
of the brackets) and the ladder operator definition of squeezed states (the curly brackets) 
p|. This completes the transformation of a coherent state in the (A, A''') representation 
into a squeezed state in the (a, a)) representation. Since the squeeze operator is invertible, 
we can obviously reverse this process. R 



6 Conclusion 

The coherent states of {Z,P) or {A,A^) are squeezed states of {z,p) or (a, a"!"), and vice 
versa. This makes sense. The fundamental potentials are of different widths, so their 



coherent-state Gaussians are also of different widths. Indeed, Eq. (|47| ) shows this. For 
the {z,p) uncertainty relation, [—ie(t) / e(t)] is a squeeze factor. 



^ In Eq. (|65|), a is indeed the quantity of Eq. (pq). However, because of our construction, it can be 
obtained from Eqs. (|6^) and (|3|) as a functional of A and /3. 

^ Because of the isomorphism of the two Heisenberg-Weyl algebras, {a,a\l} and {A{t),a^{t),I}, 
the converse of this result follows. Starting with a coherent-state equation in the (a, a^) representation: 
aDo (77)10) = r]Da{r])\0), where 77 is complex, we can map this equation into an analogous squeezed- 
state equation in the representation with an appropriate squeezing operator in the {A,A'') 
representation. This is done going through an analogous procedure as above with the aid of the equations 
a = pA + uA\ P = \ - «e*(t)) , and <J ^ \ (e(t) - ii{t)) and the condition |crp - |pp = 1. 
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